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Abstract 


The  solubility  isotherm  for  CaHP04  *2H20,  DCPD,  in  the  three- 
component  system  Ca ( 0H)2 -H3 P04 -H2 0 was  determined  in  the  pH 
range  3-5  to  6.8  at  37-5°C  by  leaching  a thermostated  column 
of  DCPD  with  dilute  phosphoric  acid  solutions.  In  confirma- 
tory experiments,  equilibrium  was  approached  both  from  super- 
and  under-saturation  by  shaking  DCPD  with  appropriate  solu- 
tions. The  calculated  ionic  activity  product,  (Ca++)  x (hpo4), 
appeared  to  be  a parabolic  function  of  pH  with  a minimum  near 
pH  5.0.  The  pH  dependence  of  the  ionic  product  could  be 
accounted  for  by  considering  the  ion  pairs  [CaHP04]°  and 
[CaH2P04]4'  as  semi-empirical  parameters.  Under  the  condition 
of  saturation  with  respect  to  DCPD,  the  activity  of  the  pair 
[CaHP04]°  must  be  a constant.  The  activity  of  the  species 
[CaH2P04]+  was  shown  to  vary  directly  with  hydrogen  ion 
activity.  The  activities  of  the  two  ion  pairs  were  adjusted 
to  give  a set  of  pH-independent  ionic  activity  products  with 
a mean  of  2.19  ± 0.011  x 10-7 . The  stability  constants  for 
[CaHP04]°  and  [CaH2P04]+  are  5.88  ± 0.031  x 102  and  7.49  ± 
0.039,  respectively.  Experiments  were  conducted  to  study  the 
hydrolysis  of  DCPD  to  more  basic  calcium  phosphates  and  the 
kinetics  of  these  transformations  is  discussed.  The  signifi- 
cance of  the  ion  pairs  in  human  serum  is  considered. 


1.0  Introduction 

Dicalcium  phosphate  dihydrate,  (DCPD),  CaHP04*2H20,  has  been  found  as  a 
component  of  dental  calculi  (l)  and  other  biological  deposits  (2).  A role 
for  DCPD  as  a nucleating  agent  (3)  in  calcification  of  tissues  and  its 
presence  at  the  dental  caries  site  (4)  have  also  been  suggested.  Also,  DCPD 
has  chemical  properties  which  are  very  helpful  for  the  study  of  chemical 
interrelations  among  the  various  phosphates  of  biological  significance;  its 
high  rate  of  solution  for  example,  makes  it  possible  to  prepare  reproducible 
calcium  phosphate  solutions  within  a wide  range  of  concentrations  and  pH 
values  which  can  be  used  in  connection  with  crystal  growth  studies,  calcifi- 
cation, etc.  It  has  commercial  importance  as  a constituent  of  fertilizers, 
feed-grade  dicalcium  phosphate,  and  dentifrices. 

For  these  reasons,  its  solubility  is  of  importance;  there  is  disagree- 
ment, however,  in  reported  values  of  its  solubility  product  constant  (5,  6, 

7,  8,  9).  Perhaps  the  most  important  cause  for  inconsistent  reports  is  its 
metastability  with  respect  to  the  anhydrous  salt,  CaHP04,  and  the  more  basic 
calcium  phosphates,  hydroxyapatite  (HA),  Ca10 (0H)2 ( P04 )s , and  octacalcium 
phosphate  (OCP),  Ca8 H2 ( P04 )6  • 5H2 0,  (10).  Formation  of  complexes  or  ion  pairs 
between  the  calcium  and  orthophosphate  ions,  as  reported  here,  also  contri- 
butes to  uncertainty  in  solubility  product  constants  calculated  from  solubili- 
bility  data.  Its  solubility  isotherm  in  the  three-component  system,  Ca(0H2 
-H3P04-H20,  at  25°C  is  known  (9).  The  solubility  isotherm  and  solubility 
product  at  physiological  temperature,  37.5°C,  are  reported  here.  Equilibra- 
tions of  DCPD  with  dilute  phosphoric  acids,  using  the  leaching  technique 
previously  described  (9),  proved  successful  in  minimizing  formation  of  secon- 
dary phases.  Batchwise  equilibrations  were  used  to  verify  the  results  ob- 
tained with  the  leaching  technique. 
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2.0  Methods  and  Experimental  Materials 

2.1  Preparation  of  Materials 

Dicalcium  phosphate  dihydrate  was  prepared  by  ammoniating  in  an  ice  bath 
an  aqueous  solution  initially  saturated  with  Ca(H2P0^3 *H20  and  CaHP04»2H20; 
details  are  described  elsewhere  (9).  Chemical  analysis  of  this  salt  yielded 
Ca,  23.20$  and  P,  17.93$  (theoretical:  Ca,  23. 20$,  and  P,  18.02$).  Petro- 

graphically,  the  preparation  appeared  as  pure,  homogeneous  DCPD. 

Phosphoric  acid  solutions  were  prepared  from  doubly  crystallized  hemi- 
hydrate,  2H3P04-H20,  and  freshly  boiled  conductance  water.  The  molarities  of 
the  final  solutions  were  checked  by  chemical  analysis. 

2.2  Leaching  Experiments 

The  leaching  apparatus  is  shown  in  Figure  1.  An  8-cm  layer  of  DCPD 
crystals  was  poured  on  top  of  the  fritted  glass  (medium  porosity),  F;  phos- 
phoric acid  solution  was  poured  into  the  column;  the  solid  was  stirred  to 
remove  air  bubbles;  finally  about  300  ml  of  phosphoric  acid  solution  was 
poured  into  the  reservoir,  S,  continuously  flushing  the  apparatus  with 
nitrogen.  Water,  from  a constant  temperature  bath,  was  pumped  through  the 
jackets,  W,  maintaining  the  temperature  at  37*5  ± 0.1°C  as  measured  by  ther- 
mometer T.  The  effluent  from  the  column  was  collected  in  the  thermostated 
cup  D in  which  were  mounted  glass  and  calomel  electrodes  of  a commercially 
available  pH  meter  (claimed  relative  accuracy  ± 0.0037  pH  units)  arid  through 
which  flowed  a slow  current  of  C02-free  and  water-saturated  nitrogen.  When 
the  pH  of  the  solution  entering  the  cup  D became  constant,  several  samples 
for  analysis  were  collected  through  stopcocks  C or  E.  The  pH  was  rechecked 
several  times  during  sampling  and  at  various  flow  rates  to  be  sure  that  the 
effluent  was  saturated.  Standardization  of  the  pH  meter  was  made  with  certi- 
fied NBS  buffer  standards  (ll).  The  uncertainty  in  pH  measurements  was 
estimated  to  be  about  ± 0.008  pH  units.  The  hydraulic  head  was  kept  constant 
by  means  of  the  Mariotte  valve  M.  The  rate  of  flow  at  which  the  samples  were 
collected  was  about  25  ml  per  hour.  When  faster  rates  were  needed,  nitrogen 
pressures  of  up  to  3 psi  were  used  to  supplement  the  hydraulic  head. 

2.3  End-Over-End  Equilibrations 

A series  of  equilibrations  was  made  in  glass-stoppered  bottles,  each  of 
which  contained  5 g of  DCPD  and  120  ml  of  phosphoric  acid  (1.00  x 10-4 , 

1.00  x 10-3,  5*00  x 10“3,  or  5.00  x 10-2M)  by  slowly  rotating  end-over-end 
in  a water  bath  at  37.5  ± 0.1°C,  In  addition  to  the  DCPD,  0.1  g of  octa- 
calcium  phosphate,  Ca8 H2 ( P04 )6  • 5H2 0,  (0CP),  was  added  to  replicates  of  the 
two  most  dilute  systems.  At  selected  time  intervals  the  contents  of  the 
bottles  were  filtered  through  the  sintered  glass  in  the  thermostated  appara- 
tus shown  in  Figure  1;  pH  values  were  determined  and  samples  were  collected 
for  calcium  and  phosphorus  analyses  in  the  manner  described  above. 

Equilibrium  was  also  approached  from  supersaturation,  taking  advantage 
of  the  negative  thermal  coefficient  of  solubility  of  DCPD:  a column  of  DCPD 

in  the  apparatus  shown  in  Figure  1 was  leached  with  5-00  x 10-3M  phosphoric 
acid  at  25  ± 0.1°C;  the  pH  of  the  effluent  was  determined  and  a sample  was 
taken  for  calcium  and  phosphorus  analyses.  The  effluent  was  collected  under 
a C02-free  atmosphere  in  a glass  bottle  containing  about  2 g of  DCPD;  when  the 
volume  was  about  120  ml,  the  bottle  was  stoppered  and  equilibrated  in  a water 
bath  at  37.5  ± 0.1°C  as  described  above.  The  pH  and  the  calcium  and  phos- 
phorus concentrations  of  the  liquid  phase  were  determined  after  16  days  of 
equilibration. 
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2.4  Calcium  and  Phosphorus  Determinations 

Calcium  was  determined  colorimetrically  by  the  indirect  technique  de- 
scribed by  Banerjee  et_  al  (12).  Phosphorus  was  determined  using  the  vanado- 
molybdate  reagent  of  BraFson  e_t  al  (13).  In  both  analyses.,  samples  were  com- 
pared directly  with  blanks  and  standards  in  a double  beam  spectrophotometer. 
Uncertainties  in  phosphorus  and  calcium  determinations  were  estimated  as 

1.0  and  2.5  percent  of  the  amounts  present,  respectively. 

3.0  Calculations 


3.1  Apparent  Solubility  Product  Constants,  K'Sp 

The  ionic  strength,  p.,  is  defined  in  these  calculations  by 

P = P[v  + H+/2P  +1/1]  [1] 

in  which  P and  H+  represent  molar  concentrations  of  phosphorus  and  hydrogen 
ion  and  the  other  symbols  are  defined  as: 

v = 2 M/P 

i * (Hi>-.  + i 

2KS  f 1 f 

X = (HP3  + ( H + ) + 1 

Ki  K2  Kaf  f" 

In  equations  [2],  [3],  and  [4],  (H+)  and  M represent  the  hydrogen  ion 
activity  and  the  molar  concentration  of  calcium,  respectively,  Kt  and  K2  are 
the  first  and  second  dissociation  constants  for  phosphoric  acid,  and  f'  and  f" 
are  the  molar  activity  coefficients  for  the  species  H2P04~  and  HP04=,  respec- 
tively. The  numerical  values  for  and  K2  used  here  were  6.53  x 10-3  (l4) 
and  6.57  x 10-8 , respectively.  The  value  of  K2  at  37-5°C  was  interpolated 
by  means  of  the  equation 

pK2  = 7-1994  - 1.135  x 10-2  X + 1.998  x 10-3  X2- 

6.64  x 10-5  X3  with  X = (T°K  - 298.l6)/5 

derived  from  a portion  (0°-45°C)  of  the  experimental  data  of  Bates  and  Acree 
(15). 

Equations  [1],  [2],  [3]  and  [4]  are  consistent  with  the  assumptions  that 
the  only  ionic  phosphate  species  present  in  significant  concentrations  are 
H2P04-  and  HP04~  and  that  the  contribution  of  hydroxyl  ion  to  the  ionic 
strength  is  negligible  in  comparison  with  that  of  the  other  ionic  species 
present.  These  assumptions  are  justifiable  in  the  pH  range  of  the  present 
investigation . 


[2] 

[3] 

[4] 


An  approximate  ionic  strength  was  obtained  by  the  use  of  equations  [l] 
to  [4]  assuming  unitary  values  for  f'  and  f",  improved  ionic  activity 
coefficients,  f^,  were  then  calculated  from  the  extension  of  Debye-Huckel 
limiting  law. 


A ZiVm. 

log  fi  = 1 + Ba-j/p 


[51 
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in  which  represents  the  valence  of  the  ith  ion,  a,  its  distance  of  closest 
approach,  and  A and  B are  temperature-dependent  constants ; the  numerical  values 
chosen  (16)  for  A and  B were  0.522  and  0.331  x 10s , respectively,  the  numeri- 
cal values  for  a^  were:  Ca++,  6 x 10-8cm;  HP04  = and  H2P0j,4  x 10“8  cm  ( 17 ) • 

With  activity  coefficients  obtained  by  the  use  of  equation  [5],  a new  value 
for  p was  calculated  from  equations  [1[  to  [4[  and  iteration  was  continued 
until  convergence  was  attained,  which  usually  occurred  after  three  cycles. 

The  activity  of  HP04~,  (HP04  = ),  was  calculated  from 

(HP04=)  = P/X  [6] 

and  a solubility  product  for  dicalcium  phosphate  dihydrate,  K^p,  was  calcu- 
lated according  to  equation  [7]*  assuming  a value  of  unity  for  the  activity 
of  water, 

K'p  - M / ca  (HP0‘=)  = M P [ca  [7] 

X 

The  solubility  product  calculated  in  this  way  was  found  to  be  dependent  on 
concentration  and  is  referred  to  here  as  the  "apparent"  solubility  product. 

A more  nearly  constant  value,  Ksp,  was  obtained  by  allowing  for  the  presence 
of  the  ion  pairs,  [CaHP04]°  and  LCaH2P04]  + as  described  below. 

3.2  Stability  Constants 

The  several  stability  and  ionization  constants  considered  in  this  paper 
are  defined  as  the  equilibrium  constants  of  the  following  reactions  in  solution: 

a)  Associative  reactions  (stability  constants) 


Ca++  + HP07  - CaHP04  [k£]  [8] 

Ca++  +■  H2P07  * Ca(H2P04)+  [K*]  [9] 

b)  Dissociative  reactions  (ionization  constants) 

Ca(H2P04)+  52  CaHP04  + H+  [Ky]  [10] 

CaHP04  = CaP07  + H+  [Kx]  [11] 

For  systems  saturated  with  respect  to  DCPD,  the  concentration  of  the 
uncharged  species  [CaHP04]°,  X°,  is  given  by 

X°  = K*  (Ca++)  (HP04  = )//xO  = K*  Ksp//xo  [12] 


in  which  (Ca++)  is  the  activity  of  ionic  calcium  and  Ksp  is  the  activity 
product  (Ca++)  (HP04=) . Therefore,  X°  is  constant,  assuming,  as  customary, 
that  the  activity  coefficient,  /x°,  of  an  uncharged  species  in  a dilute 
electrolyte  solution  is  unity.  The  activity,  (yt),  of  the  charged  species 
[CaH2P04]+  is  given  by 

(y  + ) = y+  fy+  = X°(H  + )/Ky  [13] 

in  which  Y+  represents  the  concentration  of  the  species  [CaH2P04]+;  accordingly, 
the  activity  of  the  charged  species  is  proportional  to  the  activity  of  hydro- 
gen ion.  In  the  calculations  reported  here,  it  was  assumed  that  /y  + = / . 

It  follows  then,  that  if  any  two  solutions  designated  by  "a"  and  "b"  are  con- 
sidered, the  concentration  Y-fe  in  solution  "b"  is  related  to  that  in  solution 
"a"  by  the  equation 

Yb  = (H+)b  /a  Y+ 

(H+)a  fi  3 
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For  the  calculations  in  this  paper,  solution  "a"  was  taken  as  the  most 
dilute  solution  in  Table  1.  The  selection  of  values  of  X°  and  Y+-  was  made 
from  an  array  of  reasonable  pairs  set  up  for  exploratory  purpose!,  the 
choice  of  the  values  finally  selected  was  based  upon  considerations  outlined 
in  the  discussion. 


3.3  Corrected  Solubility  Product  Constant, 

bjJ 

The  ionic  strength  was  now  calculated  considering  the  balance  equations 
[15]  to  [19] 


T = Y+  f X° 

[15] 

M = [ Cct  + ] + T 

[16] 

P = X Pu  ' + T 

J ~ 0 

[17] 

[Ca++]  = M(  l-(3 ) 

[18] 

x p(J  )=  p (1-1  pv) 
3 

[19] 

In  equations  [15]  to  [19],  t is  the  total  concentration  of  calcium  (or 
phosphorus)  in  ion  pairs,  [Ca-1"1']  is  the  concentration  of  unpaired  calcium; 
p(j)  is  the  concentration  of  unpaired  phosphate  species  of  charge  j;  (3,  the 
fraction  of  total  calcium  in  paired  forms,  is  P = t/M.  The  ionic  strength, 

\i , is  given  by 

p = P[v  - 3 (1  Pv)  + ($A)  (1-1  Pv)  - ( l/2  P ) (X°  - H+)  ] [20] 

2 2 2 

and  the  solubility  product  constant,  Kor.  is  given  by 

Ksp  = MP(l-p)  (1-1  Pv)  (fCa/x)  [21] 

2 

For  each  selected  pair  of  values  for  X°  and  Y+,  equation  [20]  was  solved 
by  successive  approximations  starting  with  unStary  activity  coefficients 
analogously  to  the  procedure  outlined  above. 

4.0  Results 

The  composition  of  effluents  in  the  leaching  experiments  and  the 
corresponding  calculated  solubility  products  are  shown  in  Table  1.  All 
the  compositions  are  averages  of  at  least  two  samples.  Within  experimental 
error,  the  final  phosphorus  concentration  equaled  the  original  phosphoric 
acid  concentration  plus  the  calcium  concentration  in  the  effluent,  indica- 
ting that  congruent  dissolution  took  place  in  each  case. 

The  calcium  and  phosphorus  concentrations  in  the  various  effluents  at 
37.5°C  are  plotted  against  pH  in  Figure  2;  corresponding  curves  for  25°C 
(9)  are  shown  as  broken  lines.  It  is  apparent  that  DCPD  is  more  soluble  at 
25 °C  than  at  37-5°C. 

Apparent  solubility  products,  K'  , listed  in  Table  1,  are  plotted 
against  pH  in  Figure  3.  The  concavespshape  of  this  curve  is  consistent 
with  a model  in  which  the  ion  pairs  [CaHP04]°  and  [CaH2P04]  + are  present 
in  solution.  The  corrected  solubility  product  constants,  Kgp,  calculated 
on  the  basis  of  this  model  are  also  plotted  in  Figure  3.  Maximum  variation 
in  KSp,  about  4 percent,  is  satisfactory  in  view  of  the  analytical  uncertain- 
ties. In  contrast,  K^p  varied  by  about  20  percent. 


5 


In  Table  2 are  listed  the  parameters  used  in  the  calculations  along  with 
pertinent  stability  and  ionization  constants.  The  values  for  K*  and  K*  were 
calculated  using  equations  [12]  and  [13],  respectively.  The  values  for  K* 
were  calculated  from  K*  = (Y+)K2 

y Ksp(H+)  [22] 

The  constancy  in  the  value  of  K is  imposed  by  the  calculative  procedure, 
as  will  be  explained  later.  The  near  constancy  in  the  values  of  K*  and  K* 
(Table  2)  lends  strong  support  for  the  existence  of  the  complexes  postu-  y 
lated  here. 

The  fraction  of  phosphorus  complexed,  t/P,  is  also  plotted  in  Figure  3. 
There  appears  to  be  a minimum  of  complexing  at  pH  5.  As  might  be  antici- 
pated , the  minimum  in  the  fraction  of  P present  as  ion  pairs  is  nearly  co- 
incident with  the  minimum  in  the  K'  curve  in  the  same  figure. 

sp 

The  compositions  of  all  the  solutions  obtained  in  the  end-over-end 
equilibrations,  listed  in  Table  3 and  shown  in  Figure  2 as  triangles,  lie 
close  to  the  isotherm  obtained  in  the  leaching  experiments.  However,  only 
two  of  the  systems,  the  last  two  entries  in  Table  3,  showed  congruent 
dissolution.  These  contained  only  DCPD  and  were  equilibrated  with  5- 00  x 
10-3M  and  5-00  x 10-SM  phosphoric  acid,  respectively.  The  compositions  of 
these  two  solutions  are  the  only  ones  that  are  essentially  the  same  as  those 
obtained  in  the  corresponding  leaching  experiments.  Table  1. 

In  the  OCP-free  systems  equilibrated  with  1.00  x 10“4M  phosphoric  acid. 
Table  3,  a more  basic  phosphate  precipitated  as  evidenced  by  the  changes 
in  composition  with  time.  The  rate  of  this  precipitation  was  relatively 
fast  at  the  beginning  and  then  decreased. 

The  two  systems  containing  both  DCPD  and  OCP,  Table  3,  were  included 
with  the  hope  of  obtaining  the  singular  point  for  the  two  salts  (i.e.,  the 
solution  which  would  be  in  equilibrium  with  both  salts).  The  lack'  of  agree- 
ment in  the  final  composition  of  the  two  solutions  indicates  that  at  least 
one  of  them  did  not  remain  at  the  singular  point  composition.  Hydrolysis 
of  OCP  to  a more  basic  phosphate,  presumably  hydroxyapatite,  would  account 
for  this  result. 

The  initial  and  final  compositions  in  the  experiment  in  which  equilibrium 
was  approached  from  supersaturation  are  shown  in  the  tabulation: 


Ca 

p 

pH 

o 

i — i 

X 

s 

M x 104 

Initial 

5.202 

52.98 

101.9 

Final 

5.142 

48.22 

100.8 

The  initial  and  final  compositions  are  shown  as  squares  in  Figure  2.  It  is 
apparent  that  the  former  lies  on  the  25  °C  isotherm  and  the  latter  on  the 
37.5  °C  isotherm. 


5.0  Discussion 
5.1  Leaching  Experiments 

The  more  basic  calcium  phosphates  (such  as  HA  and  OCP)  were  not  formed 
in  the  leaching  experiments  as  long  as  the  phosphoric  acid  solutions  were 
more  concentrated  than  1.00  x 10-4M.  Attempts  to  attain  equilibrium  with 
1.00  x 10-4M  acid  failed  because  considerable  hydrolysis  of  DCPD  was  evident 
even  at  flow  rates  3 to  4 times  higher  than  the  rate  used  with  more  concen- 
trated acids.  This  behavior  is  in  marked  contrast  to  that  at  25°C  (9)  which 
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showed  no  evidence  of  hydrolysis  with  the  same  acid. 

The  following  facts  show  that  equilibrium  was  closely  approached 
in  the  experiments  reported  here  : a)  the  plots  of  experimental  points 

yielded  smooth  curves;  b)  the  end-over-end  experiments  (systems  8 and  9, 

Table  3)*  which  represent  much  longer  equilibration  times.,  gave  solution 
compositions  in  agreement  with  those  obtained  in  the  leaching  experiments; 
c)  the  system  in  which  equilibrium  was  approached  from  supersaturation., 
agreed  with  that  in  which  equilibrium  was  attained  from  undersaturation 
with  the  same  acid;  d)  within  experimental  error,  a constant  solubility 
product  was  obtained  when  corrections  for  the  presence  of  complexes  in  so- 
lution were  introduced. 

Dicalcium  phosphate  dihydrate  is  metastable  with  respect  to  anhydrous 
dicalcium  phosphate  over  the  entire  range  studied  and  with  respect  to  OCP 
and  HA,  in  part  of  the  range.  The  material  in  the  column  was  examined 
after  each  experiment  with  the  petrographic  microscope.  With  one  exception, 
the  system  in  which  the  most  concentrated  acid  was  used,  no  anhydrous  di- 
calcium phosphate  was  found.  The  trace  found  in  the  one  exception  is  be- 
lieved to  be  too  small  to  influence  the  corresponding  point  on  the  isotherm 
because  it  is  known  that  the  rate  of  precipitation  of  the  anhydrous  salt  is 
much  less  than  the  rate  of  dissolution  of  DCPD  (l8). 

5.2  End-Over-End  Experiments 

Since  all  the  compositions  shown  in  Table  3 fall  on  the  DCPD  isotherm 
(Figure  2)  obtained  in  the  leaching  experiments,  the  rate  of  precipitation 
of  more  basic  calcium  phosphates  (such  as  OCP  and  HA)  in  the  end-over-end 
experiments  must  have  been  significantly  lower  than  the  rate  of  dissolution 
of  DCPD. 

With  only  DCPD  and  OCP  present  as  solid  phases,  the  solution  composi- 
tion should  approach  the  DCPD-OCP  singular  point.  The  formation  of  HA  in 
such  a system  would  cause  compositions  to  shift  toward  the  more  acid  region. 
This  would  create  a condition  of  undersaturation  with  respect  to  OCP;  the  OCP 
would  then  dissolve  leading  to  supersaturation  with  respect  to  DCPD.  The 
dissolution  of  OCP  and  precipitation  of  DCPD  would  thus  tend  to  maintain 
the  pH  at  a higher  level,  approaching  that  of  the  OCP-DCPD  singular  point. 
Therefore,  depending  upon  the  relative  rates  of  these  three  reactions,  the 
solution  composition  would  be  either  at  the  DCPD-OCP  singular  point  or 
somewhere  on  the  acid  side  of  this  point.  In  agreement  with  these  considera- 
tions, the  1.00  x 10-3M  acid  yielded  an  equilibrium  pH  of  6.295  in  the 
leaching  experiments;  the  same  acid  after  288  hours  in  the  end-over-end 
experiments  gave  a pH  of  5-653.  The  difference  must  be  attributed  to  the 
formation  of  HA.  Likewise,  it  is  probable  that  the  DCPD-OCP  singular  point 
has  a pH  above  the  value  5-984  obtained  with  the  1.00  x 10“4M  acid  after 
264  hours. 

The  ionic  product  for  OCP  (Ca4'"1')4  (H)  ( PO4  )3  calculated  for  the  solution 
with  pH  5-984  ((3  = 0.08)  is  3-8  x 10-48.  Because  of  the  foregoing  reasons, 
the  solubility  product  of  OCP  at  37-5°C  must  be  greater  than  3-8  x 10“48 
which  puts  it  in  the  neighborhood  of  the  product  reported  for  25°C,  1.25  x 
10"47  (10). 

In  the  absence  of  OCP,  the  1.00  x 10-4M  acid  gave  solutions  with  pH 
values  below  5-984,  the  value  obtained  in  the  presence  of  OCP;  at  24  hours 
the  pH  was  5.883  and  it  decreased  to  5.281  after  960  hours.  Precipitation 
of  HA  liberates  acid  which  is  not  neutralized  by  dissolution  of  OCP  when 
the  latter  is  absent.  Consequently,  the  pH  in  the  time  series  with 
1.00  x 10_4M  acid  in  the  absence  of  OCP,  decreased  fairly  rapidly  to 
approach  that  of  the  DCPD-HA  singular  point.  The  results  suggest  that  the 
pH  of  the  DCPD-HA  singular  point  is  below  5.281,  and  probably  above  5.132, 
the  value  obtained  in  system  8,  Table  3,  in  which  no  hydrolysis  was  detected 
after  192  hours. 
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The  hydrolysis  steps  outlined  in  the  foregoing  discussion  are  consistent 
with  the  known  behavior  of  DCPD  (9,  10)  and  of  OOP  (19).  The  formation  of 
OOP  by  hydrolysis  of  DCPD  in  aqueous  solutions  occurs  at  25°C,  but  it  was 
possible  (10)  to  determine  the  singular  point  for  these  two  salts,  because 
at  25 °C  hydrolysis  of  OOP  to  HA  is  too  slow  to  displace  the  composition 
from  that  of  the  singular  point. 


5.3  Ion-Pair  Formation 

The  values  for  the  concentrations  of  the  ion-pairs  [CaHP04]°  and 
[CaH2P04]+,  X°  and  Y+,  in  the  most  dilute  solution  were  selected  from  an 
array  according  to  the  following  criteria: 


a)  the  calculated  value  for  K cannot  be  higher  than  the  minimum 
value  for  K'  -in  Table  1;  b)  the  sum  of  the  concentrations  of  the  ion  pairs 
cannot  be  higher  than  either  the  total  Ca  or  P concentrations  in  any  of 
the  systems  considered;  c)  there  must  be  a minimum  in  the  fraction  (t/p)  of 
P complexed  and  it  should  coincide  with  the  minimum  in  the  value  of 
in  Figure  3 ; d)  the  variance  of  Ks„  obtained  with  the  final  values  of  X° 
and  Yg  should  be  an  absolute  minimum  with  respect  to  all  other  possible 
pairs  of  values  for  these  parameters. 

Designating  values  for  X°  and  Y’1’  for  the  most  dilute  solution  fixes 

a numberical  value  of  K for  all  the  other  solutions  because  the  quantities 

(Y+)/(H+)  and  X°  in  equation  [13]  are  constants  in  the  adopted  model.  This 

is  why  the  quantity  Ky  does  not  appear  in  Table  2 along  with  the  calculated 

values  for  K*  and  K*.  The  value  derived  for  Klr  is  5'.  16  x 10-6  . 

j y 

The  value  for  the  acid  ionization  constant,  K , of  the  species 
[CaH2P04]^,  5.16  x 10  s,  is  one  hundred  times  the  Second  ionization  constant 
of  phosphoric  acid;  this  ion-pair  is  thus  a stronger  acid  than  H2P04,  pre- 
sumably due  to  the  proximity  of  the  charge  of  the  calcium  ion. 

If  the  reasonable  assumption  is  made  that  the  ratio  of  IC./K  is  of  the 
same  order  of  magnitude  as  Kj/Kp  or  K2/K3  for  phosphoric  acid,  the  value  of 
Kx  would  be  of  the  order  of  10~ri . This  means  that  the  activity  of  the 
species  [CaP04]“  in  the  system  with  the  highest  pH  would  be  of  the  order  of 
10-SM,  which  is  neglibible  relative  to  the  activities  of  Ca+(J[  and  HP04“. 
Therefore,  ignoring  the  concentration  of  the  species  [CaP04]-  in  the 
calculations  is  justified. 

In  a previous  study  of  the  solubility  of  DCPD  at  25°C  (9)j  the 
possibility  of  complexing  was  not  ruled  out,  but  the  data  were  treated  by 
an  empirical  method  that  gave  K by  extrapolation  to  zero  ionic  strength. 
The  25°  data  are  being  re-evaluated  in  this  laboratory,  taking  into  account 
ion-pair  formation. 


According  to  the  present  treatment,  it  appears. that  as  much  as  20  per- 
cent of  the  calcium  was  ion-paired  (e.g.,  system 'no.  8,  Table  l).  The  model 
adopted  for  the  treatment  of  the  ion  pairs  is  not  exempt  from  uncertainties; 
for  example,  it  is  impossible  to  ascertain  the  adequacy  of  equation  [5]  for 
calculating  activity  coefficients.  On  the  other  hand,  the  variation  of  the 
solubility  product  constant  could  not  be  explained  simply  in  terms  of 
activity  coefficients  because  a)  the  extended  Debye-Huckel  limiting  law 
should  give  reasonably  good  values  for  activity  coefficients  in  solutions  as 
dilute  as  those  used  here,  and  b)  it  is  unlikely  that  the  minimum  in  the 
value  of  K'  , which  occurs  at  |i  « .025  is  related  to  minima  in  mean 
activity  coefficients  which  are  sometimes  observed;  such  minima  occur  at  con- 
centrations much  higher  than  those  described  here. 
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5.4  Comparison  with  Reported  Data 

Davies  and  Hoyle  (20)  reported  dissociation  constants  at  25°C  for 
systems  containing  sodium  and  potassium  phosphate  buffers  and  calcium  iodate. 
Allowing  for  experimental  differences,  the  constants  obtained  by  Davies  and 
Hoyle,  K*  (25°)  = 1/.002  = 5 x 10s  and  K*  (25°)  = 1/.084  = 11.8,  are  in 
reasonable  agreement  with  those  given  here,  K*  (37.5°)  = 5.88  x 102  and  K* 
(37.5°)  = 7.49.  X y 

Bjerrum  (21)  studied  the  three  component  system  Ca(0H)2 -H3 04 -H20  at 
l8°C  in  the  pH  range  6 to  8 in  the  absence  of  solid  phases.  He  concluded 
that  a model  involving  the  ion-pair  [CaHP04]°  or  any  of  the  ion  pairs 
[Ca(HP04  )2  ]2“,  ( Ca  ( HP04  ) (P04)]3-,  or  [Ca(  P04  )2  ]4  ",  would  fit  his  data. 

Models  involving  [CaHP04]°  alone  or  in  combination  with  any  of  the  above 
ion  pairs  with  Ca/P  ratios  of  1/2  gave  variable  Kg  values  when  used  with 
the  present  data.  The  only  model  giving  reasonably  constant  K values 
was  the  one  with  the  ion  pairs  [CaHP04]°  and  [CaH2P04]  + . It  is^believed 
that  the  condition  of  saturation  with  respect  to  DCPD,  which  demands  a con- 
stant activity  for  the  species  [CaHP04]°,  is  more  determinative  than  the 
conditions  used  by  Bjerrum. 

More  recently,  Greenwald  (22)  used  Bjerrum' s data  to  calculate  values 
for  the  instability  constant  of  [CaHP04]°,  these  varied  from  1.01  to 
9.55  x 10-3  . As  in  previous  reports  by  the  same  author  (23.,  24),  ionic 
concentrations  were  used  instead  of  activities.  For  this  reason,  and  be- 
cause the  procedure  used  in  his  calculations  is  not  detailed,  a direct 
comparison  of  his  results  with  those  given  here  is  not  possible.  His 
average  value  for  the  instability  constant,  2.77  x 10“3,  gives  a K*  (l8°C)  = 
1/2.77  x 10-3  = 3.6l  x 102 . This,  too,  gives  order  of  magnitude  agreement 
with  the  value  of  K£  (37.5°C)  = 5.88  x 102  given  here. 

Similar  rough  agreement  is  found  with  the  result  of  Gosselin  and 
Coghlan  (25),  K*  (25°C)  = I/I.38  x 10“2  = 0.72  x 102 , who  also  used  ionic 
concentrations  instead  of  activities  in  their  calculations. 

The  negative  temperature  coefficient  for  the  solubility  of  DCPD  relates 
to  properties  of  both  the  solution  and  the  solid.  For  example,  with  an  in- 
crease in  temperature  the  activity  of  water  increases  and  it  is  conceivable 
that  less  water  is  associated  with  the  calcium  and  phosphate  ions  which, 
then,  for  comparable  concentrations,  would  display  a higher  activity;  this 
should  result  in  lower  solubility  of  DCPD  as  the  temperature  increases. 
However,  if  the  solublity  product,  Ks  , decreases  as  the  temperature  in- 
creases, the  magnitude  of  the  standard  free  energy  of  formation  of  DCPD 
must  increase  with  temperature  more_rapidly  than  the  sum  of  the  free 
energies  of  formation  of  Ca++,  HP04“  and  2H20.  The  solubility  product  of 
DCPD  at  25°C  (9),  2.77  x 10~7 , is  higher  than  the  value,  2.19  x 10-7 , 
at  37*5°C.  However,  since  no  consideration  was  given  to  the  formation  of 
ion  pairs  in  solution  at  25°C,  a meaningful  comparison  is  not  possible  at 
present . 

An  important  application  of  the  stability  constants  relates  to  the 
formation  of  ion  pairs  in  biological  fluids.  Representative  figures  cited 
for  ionic  Ca++  and  HP04=  activities  in  adult  human  serum,  4.7  x 10  4 and 
1.9  x 10_4M,  respectively  (26),  give  the  result  at  pH  = 7-4  that  about 
five  percent  of  the  calcium  or  inorganic  phosphorus  is  in  the  form  of  the 
ion  pair_  [CaHP04  ] 0 . In  infant  serum,  with  Ca'1"1'  and  HP04  activities  of 
4.7  x 10  4 , the  percentage  of  calcium  and  phosphorus  in  paired  forms  would 
amount  to  8 and  5j  respectively.  For  both  adults  and  infants,  the  con- 
centrations of  the  ion  pairs  [CaH2P04]+  and  [CaP04]_  in  serum  appear  to  be 
negligible . 
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Solubility  of  CAHP04  *2H20 
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Averages  given  with  the  standard  error  in  the  mean.  The  value  for 
Kv  was  5-160  x 10-6  . 


End-Over-End  Equilibrations 


Figure  1.  Leaching  apparatus. 
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Figure  2.  Solubility  isotherm  of  DCPD  at  37-5°C,  full  line,  and 
at  25°C,  broken  line  (9)-  Open  circles,  compositions 
obtained  in  leaching  apparatus;  triangles,  compositions 
obtained  in  end-over-end  experiments;  squares,  initial 
and  final  compositions  in  super-saturation  experiments 
(see  tabulation  in  text). 
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Figure  3.  Solubility  product  constants  (apparent  and  corrected) 
and  percent  P in  ion  pairs  plotted  against  pH. 
Calculated  errors,  based  on  analytical  uncertainties, 
are  shown  for  each  point . 
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